Theory of the magnetoeletric effect in a hghtly doped high-Tc cuprate 
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In a recent study Viskadourakis et al. discovered that extremely underdoped La2Cu04+cc is a 
relaxor ferroelectric and a magnetoelectric material at low temperatures. It is further observed that 
the magnetoelectric response is anisotropic for different directions of electric polarization and applied 
magnetic field. By constructing an appropriate Landau theory, we show that a bi-quadratic mag- 
netoelectric coupling can explain the experimentally observed polarization dependence on magnetic 
field. This coupling leads to several novel low-temperature effects including a feedback enhancement 
of the magnetization below the ferroelectric transition, and a predicted magnetocapacitive effect. 

PACS numbers: 64.70.P-, 74.72.Cj, 77.80.-e, 77.80.Jk 



The field of magnetoelectrics has witnessed intense the- 
oretical and experimental progress in the recent years, 
mainly driven by the discovery of new materials exhibit- 
ing a non-linear magnetoelectric effect i.e. the dominant 
magnetoelectric coupling is of higher order than a bilin- 
ear coupling between electric and magnetic fields [l|, Q . 
Among the materials discovered with such unusual phys- 
ical properties are the so-called birelaxors Q . These sys- 
tems show both relaxor ferroelectric and relaxor magnetic 
properties and are associated with spin-charge coupling 
at a mesoscopic scale. Focussing on the parent high- 
Tc superconductor La2Cu04+aj (LCO) with exception- 
ally low carrier concentration n = lO^^cm"^, we have 
recently found that this material is in fact a ferroelec- 
tric at low temperatures [J]. More specifically, LCO has 
been shown to be a relaxor ferroelectric where the dielec- 
tric mode behavior is caused by freezing of randomly ori- 
ented polarized regions [5|,[6[. In addition, LCO exhibits a 
distinct magnetoelectric effect with a pronounced depen- 
dence of the polarization on an externally applied mag- 
netic field. These results shine new light on the nature of 
impurities and doped charge carriers in antiferromagnetic 
Mott insulators. Here, we focus on the magnetoelectric 
effect discovered in Ref.|3, and show how non-linear terms 
coupling polarization and magnetization naturally lead 
to the observed field effect. 

Parent cuprate superconductors are 2D antiferromag- 
nets with weak interplanar exchange coupling giving rise 
to 3D long-range Neel order [7, 8]. In the case of LCO, 
the Cu spins are slightly canted out of the Cu02 planes 
because of a finite Dzyaloshinskii-Moriya (DM) interac- 
tion existing in the low temperature orthorhombic phase 
(LTO) (sf. Though this allows a small ferromagnetic mo- 
ment to build up on each Cu02 plane, the net magnetic 
moment is zero since the moments in consecutive lay- 
ers are oriented in opposite directions. On application 
of an external magnetic field a first order spin flop tran- 



sition is observed at a critical magnetic field i/gf ~ 5T 
[sl, |lo| . Clear evidence for coupled spin and charge de- 
grees of freedom in these systems come from the obser- 
vation of pronounced discontinuities in resistivity and di- 
electric constant at a magnetic field corresponding to i^gf 
0, [HI- Further evidence of such coupling is found in 
LCO by the possibility to detwin these crystals e.g. by 
the application of an in-plane magnetic field [l2| . 

In the following we use a symmetry based analysis to 
identify the particular magnetoelectric interaction terms 
that are responsible for spin-charge coupling in under- 
doped LCO. We further construct a Landau theory, and 
show how this model reproduces all the qualitative fea- 
tures of the magnetic field dependence of the polarization 
curves reported in Ref. 0. 

Below approximately 530K the crystal structure of 
LCO is LTO with space group Cmca (D2^). If {x^y^z) 
denote unit vectors along the crystal axes then we define 
fractional translations by 



r = i(ax + cz)r = ^{ax + by), 



(1) 



where a, 6, c are the lattice constants. The symmetry 
elements of this crystal structure are then written as G = 
Go + tGo, where Go contains the eight elements 

/, da, (j'l), (j'c, C2a, ^^25, ^2^ (2) 

Here E denotes the identity, I inversion about a Cu 
site, cra,cr5,(Jc reflections about the planes x = a/2,y = 
b/2,z = c/2, and C2a-,G2h-,G2c are 180^ rotations about 
the axes that emanate from the center of the unit cell. 
Primed elements must be complemented by translation 
t' that is itself not a symmetry operation. 

Taking into account the above mentioned symmetry 
properties of the Cmca space group, the free energy can 
be expressed as a sum of three contributions 



F = Fm + Fmp + Fp . 



(3) 
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Here, Fm is the purely magnetic free energy, Fmp is the 
magnetoelectric contribution, and Fp is the polarization 
free energy. The magnetic free energy that accounts for 
the crystal structure of the LTO phase has been studied 
previously by e.g. Thio et al. [13] and is given by 



FM = \jZV^L^ + \ALt + \BLt-CL 



Xo 



, ^ Mf-H,M,-HabL,]^-J^L,L2. (4) 

Here, the out-of-plane (c direction) [in-plane {a—b plane)] 
applied magnetic field is represented by He [Hab]- The 
order parameter Mi = {SAi-^-SBi)/"^ is the ferromagnetic 
moment per spin with Sai^ Sbi being the sublattice spins 
in the i^^ plane, and Li = {Sai — Ssi)/"^ is the antiferro- 
magnetic order parameter (L^||a). The spins are slightly 
canted due to the DM interaction term — CM^L^, which 
causes them to lie in the a — c plane of the magnetic 
unit cell. The coupling between the different planes is 
included by the J± term. 

The presence of an inversion symmetry in the space 
group of the crystal forbids any linear magnetoelectric 
effect [1^ and the physics is dominated by non-linear 
coupling terms. We can focus on the largest non-linear 
terms by further noting that the experimentally observed 
polarization response is symmetric under inversion of the 
external magnetic field (i.e P{H) = P{—H)). This im- 
plies that the dominant couplings are of even order in the 
magnetic order parameter. Hence, the following terms 
contribute to the magnetoelectric coupling 



Fmp 



72a 



+73«M,L,)Pj, (5) 



where the components for P run over a = (a, 6, c) in the 
magnetic unit cell. The jia and 720; terms have been 
introduced using symmetry arguments alone but their 
microscopic origin can e.g. originate from 



Hint — —^me SiSjak^'i. 



(6) 



ij kl 



This Hamiltonian describes a bi-quadratic coupling be- 



tween spins Si J and structural pseudospins ak,i P^. A 
bi-quadratic coupling term has also been derived by Pirc 
et al. [l6|. In multi-glass material like doped SrTiOs 
fl7| or in EuTiOs [l8| these terms have been invoked to 
explain the observed magnetoelectric effect. Finally, the 
DM induced bi-quadratic magnetoelectric coupling term 
with coefficient j^a has been used to explain magneto- 
electricity in BaMnF4 

The polarization free energy is given by 



-1 

ea p 

2 ' 



13 



a 



Pa) - EP. 



(7) 



Here, Xea is the electric susceptibility for the a compo- 
nent of the polarization. Fourth order terms have been 
included to obtain stable solutions, and E denotes the 
applied electric field. In its most general form, the free 
energy should also contain gradient terms since we are 
dealing with a relaxor system [l6| as well as higher order 
terms. However, since we are interested in the magne- 
toelectric effect close to the ferroelectric transition, it is 
reasonable to restrict our analysis to the above terms. 

The solutions that determine P{H) are obtained by 
minimizing F with respect to the electric polarization 
and magnetic order parameters. In the case of LCO stud- 
ied experimentally, the Neel temperature is Tn ~ 320K, 
which is much higher than the temperature at which the 
ferroelectric order sets in (Tp ~ 4.5K) Therefore, we 
evaluate Fm for the high temperature phase with P = 0, 
providing the following set of equations 



Mi = xoiHc + CLi), 



(8) 



X2hLi + AL\ + BL\ + -J^L2 = CM^+Hab, (9) 

X2dL2 + All + +BLl + i J^Li = CM2 + Hab, (10) 
2 

IXea + Y.^llaL} + 72aM2 ^ ^g^M^L^JPa = -pPl 
i=l 

(11) 

The experimental magnetization curves at low tem- 
peratures (<~ 30K) have a glassy contribution that can 
be observed in Fig. [TJa) as a hump feature at low mag- 
netic fields and a broadened spin flop transition for out- 
of-plane magnetic fields He [20| . Such glassy behavior is 
absent in the magnetization response for in-plane applied 
magnetic fields Hah as seen from Fig. [IJb). These fea- 
tures cannot be obtained from the above equations and 
hence to include them to the lowest order, we simply take 
the experimental magnetization curves in Fig.[T]as input. 
This approximation along with the temperature range at 
which we evaluate the magnetoelectric effect leads to de- 
viations from the parameters of the magnetic free energy 
from those used e.g. by Thio et al.\l'A. |21| 

In terms of the following rescaled quantities = 
XoC(Li + L2)/2, l_ = xoC'iL, - L2)/2, M = {M^ + 





FIG. 1: (Color online) (a) Out-of-plane measured magnetiza- 
tion, and (b) in-plane magnetization measured at T = 5K. 
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^2)/2, = 27i,(xoC)-2, = 272c., 7L = 

273a ixoC) ^ , the polarization dependence on the apphed 
magnetic field can be expressed as 



Pa{0) 
Pa{Hc) 



[1 



E{0) 



{ll{Ha,)-ll{S))^ll{Ha,))]'^\ 

(12) 



[l + ^G-(^c)-/l(0)+^.M(i7,)^ 



Pa(0) 

+ (1 - - qa)ll{Hc) + qo.M[Hc)l+{HcW\ (13) 



where 5a = Ac/?_(0)/(x 
= Isa/^a with A, 



' Ac/^(0)), ^« = 7L/^«, and 
+ + tL- In general, 
all three parameters Sq,, and will be temperature 
dependent. The temperature dependence of primar- 
ily results from its relation to the electric susceptibility; 
hence an estimation of the amplitude of can reveal the 
magnitude of anisotropy in the electric polarization. 

In the case of an out-of-plane magnetic field He, the 
measured Pc{Hc) increases with field and exhibits a pro- 
nounced hump at the spin-flop transition at Hsf as seen 
from Fig. [2fa). By contrast, Pab{Hc) decreases with in- 
creasing He but also exhibits a hump feature at i^gf as 
seen from Fig. [2fc). This behavior indicates a scenario 
where the coupling between the magnetic order and out- 
of-plane polarization Pc is attractive whereas the cou- 
pling with Pab is repulsive. However, we find that the 
theoretical picture is more complex due to the presence of 
two competing magnetic orders L and M coupling to the 
electric polarization. For the out-of-plane magnetic field, 
we calculate the Pa{Hc) response using the magnetiza- 
tion data shown in Fig.flja). As seen from Fig.[2fb,d), we 
obtain qualitative agreement with both experimental po- 
larization curves including the hump feature at i^gf . Note 
that for the coupling to the magnetization we have set 
9a = 9b = 9c = —0.2. Therefore, the source of anisotropy 
between Pc{Hc) and Pab{Hc) is the DM induced coupling 
term qo^P^MiLi with = qi^ = {) and qc = —6.58. 

In the case of an in-plane magnetic field Hab, Pa{Hab) 
depends on only a single fitting parameter that con- 
trols the magnitude of the polarization ratio, whereas 
the shape of the theoretical curves are governed by the 
magnetic order parameter of the system. Further, we 
can deduce from Eq. ([12]) that the polarization primar- 
ily couples to the 2D antiferromagnetic order through a 
magnetoelectric interaction term Fmp ~ XaP^L^. The 
anisotropy in polarization through the parameter A^ for 
in-plane magnetic fields is primarily controlled by the DM 
induced magnetoelectric couphng term g^c^P^M^L^, since 
qoc undergoes the largest variation between in-plane and 
out-of-plane directions. The DM physics therefore plays 
an important role in generating an anisotropy between 
Pc{H) and Pab{H) in La2Cu04+:c. 

Using the magnetization curve from Fig. [Hb), the so- 
lutions of Pa{Hab) are presented in Fig. [3l^b,d). Though 
the theoretical curves have similar shapes, the scales are 
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FIG. 2: (Color online) Polarization Pa vs. out-of-plane 
magnetic field He. (a,c) Experimentally measured values at 
T = 2.5K. (b,d) Theoretically calculated Pa{Hc)/Pc{0) for 

Sc = 0.1, Sa = 0.074, ga = Qb = Qc = -0.2, Qc = -6.58, 

and Qa — Qb = using available experimental magnetization 
values at T = 5K. 



different in the two figures due to the difference in Sa val- 
ues. Note that all theoretical curves have been plotted 
using magnetization data available at T=5K while the ex- 
perimental curves correspond to T=2.5K which accounts 
for the smaller scale of the theoretical curves. As seen 
from the plots, we find again qualitative agreement with 




FIG. 3: (Color online) Polarization Pa vs. in-plane magnetic 
field Hab- (a,c) Experimentally measured values at T = 2.5K. 
(b,d) Theoretically calculated Pa{Hab)/ Pa{0) for Sc = 0.1 
and Sa — 0.074 using available experimental magnetization 
values at T = 5K. 
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the experimental data shown in Figs. [3ja,c) though the 
experimental plot inlSfa) has a steeper slope than theory. 
As in Fig. [21 the plots correspond to Sc = 0.1, Sa = 0.074 
which leads to a small anisotropy in the electric suscepti- 
bilities Xea therefore between the zero magnetic field 
values of in-plane polarization Pab and the out-of-plane 
component Pq. This small anisotropy has been observed 
in experiments Q and implies a much weaker anisotropy 
of the electric polarization compared to the magnetic or- 
der. This result lends support to a scenario where the 
non-stoichiometric oxygen dopants (charge carrier dop- 
ing) play an important role in generating the relaxor fer- 
roelectricity in La2Cu04+a^. 

We have observed experimentally that the magnetiza- 
tion shows a small enhancement below the temperatures 
where the ferroelectric order sets in [22| . This effect is in 
addition to the typical upturn in magnetization near the 
spin glass freezing temperature [H, ji^. We can study 
such feedback effect of a finite polarization on the magne- 
tization by minimizing F within the ferroelectric phase. 
This gives to lowest order 

XoHe + [l-XO Ea 7LPaiHc)]l+m 



Mr, 



Mab 



(14) 



(15) 



Note that in this expression the relative sign of the coef- 
ficients can be determined from the relations ")'za.h'2a. — 
Qcx/goc > 0. This also implies that since < in our 
model, it naturally causes an enhancement of magneti- 
zation due to the presence of a ferroelectric state. Ad- 
ditionally, we also find that this enhancement is present 
both in the in-plane and out-of-plane magnetization with 
the relative size of the enhancement depending on the 
amplitude of polarization change with magnetic field. 

A bi-quadratic coupling has been used to explain the 
observation of niagnetocapacitive effects in materials like 
doped SrTiOs [18|. It is defined by the relation = 
—d'^F/dP^ and hence requires at least quadratic terms 
in the polarization. For LCO the relative change in di- 
electric constant is given by 

Ae« = [cc^iH) - ec.(0)]/ec.(0) = AP^{H). (16) 

A weak magnetocapacitive effect is therefore predicted 
in experiments at low temperatures. Note that in the 
above expression we would expect a small suppression in 
the permittivity for magnetic fields in the a — b plane. 

The observation of relaxor ferroelectricity in under- 
doped LCO has been argued to originate from the for- 
mation of polar nano-regions (PNR) around the non- 
stoichiometric oxygen dopants The formation of PNR 
and the mechanism by which they condense into a ferro- 
electric phase is a well studied topic [gI, [H, [26| . Though 
the relaxor physics in LCO naturally relates to the pres- 
ence of dopants, the extremely low concentration of ex- 
cess oxygen in the samples used in Ref. .^j may imply 
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FIG. 4: (Color online) Illustration of polar nano-regions 
with randomly oriented electric polarizations (black arrows) 
at high T within the antiferromagnet (red arrows). Also 
shown are the distorted magnetic moments within a corre- 
lation length of the polar nano-regions. 



the presence of additional mechanisms for the PNR to 
couple and undergo a spontaneous transition to long- 
range ferroelectric order. One may speculate that such 
mechanisms include subtle non-centrosymmetric struc- 
tural distortions in the host lattice [lo| and/or a tendency 
for the dopants to cluster and thereby reduce the inter- 
PNR distance. As shown in Fig. IH in magnetic ma- 
terials like LCO the PNRs may also cause a distorted 
spin structure that could lead to a magnetoelectric ef- 
fect through e.g. geometric frustrations in the presence 
of DM interaction and/or indirectly through coupling to 
strains. This physics has similarities to the observation 
of magnetoelectric behavior in a number of other relaxor 
ferroelectricsllTL 1271. 



In summary, we have shown that the magnetoelectric 
effect in extremely underdoped La2Cu04+aj can be ex- 
plained by bi-quadratic terms in the free energy. It is 
proposed that the microscopic origin of the ferroelectric- 
ity is caused by polar nano-regions generated by dopant 
ions. The discovery of ferroelectricity and magnetoelec- 
tric effect in the cuprate materials due to charge carrier 
doping have sparked many questions for future studies. 
In particular, what happens at higher doping levels and 
what is the fate of the PNRs in the regime where the 
pseudo-gap and superconducting phases emerge? 
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